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Abstract: 

The microscopic quasiparticle-phonon model, MQPM, is used to study the energy 
spectra of the odd Z = 53 — 63, A = 82 isotones. The results are compared 
with experimental data, with the extreme quasiparticle-phonon limit and with 
the results of an unrestricted 2sld0gj/2^hii/2 shell model (SM) calculation. The 
interaction used in these calculations is a realistic two-body G-matrix interaction 
derived from modern meson-exchange potential models for the nucleon-nucleon 
interaction. For the shell model all the two-body matrix elements are renormalized 
by the Q-box method whereas for the MQPM the effective interaction is defined 
by the G-matrix. 
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1 Introduction 

The quasiparticle-phonon coupling (QPM) scheme is a convenient way of ex- 
ploiting the rich data on low-energy excitations of doubly-even nuclei in studying 
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nuclear-structure effects in odd-mass (odd-A for short) nuclei. Since it was intro- 
duced [m it has been used extensively, in various forms to discuss energy- 
spectra of odd-A nuclei. This scheme makes use of the property of the BCS 
quasiparticles being the elementary excitations of an odd-A nucleus and assumes 
that by coupling them to the few lowest-energy (collective) excitations of the 
doubly-even reference nucleus one is able to describe, at least qualitatively, the 
spectroscopy of odd- A nuclei. In many cases only the first quadrupole and/or 
octupole phonon have been used in the calculations. 

In the traditional quasiparticle-phonon models (QPM) independent model hamil- 
tonians have been used to create the quasiparticles, the phonons and their cou- 
pling. This means that the two-body interaction matrix elements of the many- 
fermion hamiltonian are not internally consistent. For the phonons one usually 
adopts very simple model hamiltonians, like pairing-plus-quadrupole or semiem- 
pirical ones while the quasiparticle energies can be extracted from simple pairing 
hamiltonians or directly from experiments. The quasiparticle-phonon coupling 
term has usually the simplest possible phenomenological form and may have 
adjustable parameters to control the energy of the quasiparticle-phonon multi- 
plets. In Ref. |Q a model was introduced, the so-called microscopic quasiparticle- 
phonon model (MQPM), which bears resemblance to the traditional QPM, but 
which uses a microscopic hamiltonian and a scheme to optimize the size of the 
quasiparticle-phonon basis used in the calculations. In this sense, see the discus- 
sion in section |^, it represents also an improvement over the traditional QPM 
since the quasiparticle-phonon interaction is treated in a more consistent way. 

In order to test the MQPM method, we have singled out the N = 82 isotones. The 
sequence of semi-magic = 82 nuclei shows a high degree of regularity which 
makes them well suited for systematic studies and for testing of microscopic 
nuclear models. Pairing effects seem to play an important role in the even nuclei, 
and the low-lying states in the odd nuclei may be described as one-quasiparticle 
states with increasing fermi level. Thus, models like the MQPM may be viewed 
as a reasonable starting point for the description of such nuclei. 

The first aim of this work is therefore to compare the MQPM with the more 
traditional QPM and QRPA approaches in order to see how well the experimental 
spectra can be reproduced and to interpret eventual differences. 

A comprehensive study of the even = 82 isotones was carried out by us in 
a previous work 0. There a comparison was made between the QRPA and the 
results obtained with an extensive shell model calculation. Energy spectra were 
generally well described and so were the transition probabilities. In Ref. one 
of the aims was to calculate an effective interaction based on modern meson- 
exchange models for the nucleon-nucleon (NN) potential. The first step in the 
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derivation of an effective interaction Vefr was to renormalize the NN potential 
through the so-called G-matrix. The G-matrix was in turn used in a perturbative 



many-body scheme, see e.g. Ref. |T^ for further details, to derive an effective 
interaction appropriate for the = 82 isotones. The effective interaction is meant 
to take into account degrees of freedom not included in the model space. This 
interaction was, in turn, applied in a full shell model calculation with a model 
space consisting of the orbitals 2si/2, 1(^5/25 1*^3/25 ^Qi/i and 0/;,ii/2 for the Z = 
52 - 64, = 82 isotones. 

The second aim of this work is therefore to extend the comparative analysis of 
Ref. 1^ for the even isotones to the case of the odd isotones. We will employ the 
same effective interaction and model space in the shell model analysis. Moreover, 
the same G-matrix used in the QRPA studies of Ref. ^ will be used in the 
MQPM calculations. Since the MQPM discussed here (and the QRPA method 
of as well) employs a larger single-particle space than the perturbative many- 
body scheme, our hope is to see whether the two approaches could shed light on 
different many-body contributions and their influence on various spectroscopic 
observables. 

The A^ = 82 isotones have previously been studied extensively by Heyde and 
Waroquier [pJ] , p!2| by using the quasiparticle Tamm-Dancoff approximation in- 



cluding one- and three-quasiparticle states as basis states and projecting out 
the spurious three-quasiparticle components. The calculations were done using 



the surface-delta interaction [|ri|] and a more elaborate phenomenological interac- 
tion |T2[ in the proton 2si/2lc?5/2lc?3/20(77/20/iii/2 single-particle basis, i.e. in the 
same basis which we adopt for the shell model in our present article. However, the 
more realistic results of Ref. |]T2| are not accessible for direct comparison with the 
present results of the MQPM for the following reasons: A) the MQPM employs 
a larger valence space where also neutron degrees of freedom are included. B) In 
the MQPM we use the same proton single-particle energies, extracted from the 
experimental spectrum of ^^^Sb, for all the odd isotones whereas in Ref. [O] the 



inverse-gap-equation method was used to extract single-particle energies from the 



experimental spectra of each odd N = 82 isotone separately. C) In Ref. the 
phenomenological force, namely a central force of Gaussian shape with spin ex- 
change, was fitted from case to case by the above-mentioned inverse-gap-equation 
method, and overall by level schemes and transition rates of the even A^ = 82 
isotones. Contrary to this, in the present MQPM calculation the same bare G- 
matrix interaction is used for all the even and odd N = 82 isotones without any 
fitting procedures. For the above reasons we refrain from direct comparison with 
the results of Ref. in this article and, instead, concentrate on comparison 
with the shell model results. 

The theoretical framework of the MQPM is presented in section 0. A brief review 
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of the effective interaction theory and the shell model follows in section ^. The 
results are presented and discussed in section ^. Finally, in section || we draw our 
conclusions. 



2 Theoretical framework of the MQPM 



The microscopic quasiparticle-phonon model, MQPM, represents a scheme to 
treat all the three parts of the hamiltonian, namely the quasiparticle, phonon 
and quasiparticle-phonon terms, on equal footing. This is possible by starting 
from a microscopic hamiltonian with two-body matrix elements derived from ef- 
fective matrix elements such as a G-matrix. The G-matrix is a medium modified 
nucleon-nucleon interaction where all ladder type diagrams are summed to in- 
finite order. This method enables one in a systematic way to derive both the 
proton-neutron and the like-nucleon two-body interaction. In the quasiparticle 
language these parts of the interaction relate to the H^i and H22 parts of the 
quasiparticle representation of the nuclear hamiltonian. The H22 part is treated 
in the BCS and quasiparticle random-phase approximation (QRPA) framework 
and leads to definition of the quasiparticles and the excitation (phonon) spectrum 
of the doubly-even reference nucleus. The H^i part is then diagonalized in the 
quasiparticle-phonon basis discussed below. 



The MQPM treats the structure of the odd- A nuclei in four steps. First, the 
neighboring even-even nucleus, or nuclei, can be used to study the properties 
of the chosen mass region and to fix the possible free parameters of the model 
hamiltonian. In the present C3iSG, clS also in 0, we have used a G- matrix derived 
from modern meson-exchange potential models, and thus no phenomenological 
renormalization of the two-body interacton was done. This hamiltonian is used 
to generate the phonons which are excitations of the even-even nuclei. In the 
MQPM the phonons are derived by the use of the quasiparticle random-phase 
approximation (QRPA) procedure |T3|. Second, the monopole part of the same 



hamiltonian is used to generate the quasiparticles, which are the basic building 
blocks of the odd-A excitations, through the BCS procedure. As the third step, 
the two basic excitations, QRPA phonons and BCS quasiparticles, are coupled to 
form a wave function basis for a realistic treatment of the odd-A nucleus. As the 
last step, the residual hamiltonian, containing the interaction of the odd nucleon 
with the even-even reference nucleus (the H^i + H13 part of the hamiltonian in 
Eq. (3) below) is diagonahzed in this (over-complete) basis. 
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In the MQPM the starting point is the A-fermion hamiltonian 

H = Y1 ^»clCa + jYl "^afS^/SclclcsC^, (1) 
a ^ a/375 

containing antisymmetrized two-body matrix elements Vap-yS = {aP\v\'yS) — {aP\v\S'j) 
obtained from the Bonn-A G-matrix. Greek indices denote all single-particle 
quantum numbers a = {a,ma}, and roman indices, when used, denote all single- 
particle quantum numbers except the magnetic ones, i.e. a = {na,la,ja}- 

The approximate ground state of the even-even reference nucleus is obtained from 
a BCS calculation, where quasiparticle energies and occupation factors Ua and Va 
are obtained from the Bogoliubov-Valatin transformation to quasiparticles 



= '^t^cl + v^Cf, , (2) 

where aj^ = aL^(— 1)-'^"^ and = cL^(— l)-'^"*. After this transformation the 
hamiltonian can be written in the form 

H = ^ Eattl^aa + H22 + H40 + Ho4 + Hsi + ifi3, (3) 

a 

where Ea are the quasiparticle energies and other terms of the hamiltonian are 
normal-ordered parts of the residual interaction labeled according to the number 



of quasiparticle creation and annihilation operators which they contain . 



In the first version of the MQPM ||^ the coupling part of the microscopic hamil- 
tonian, does not emerge from the equations-of- motion method (EOM) [ p!5 



The EOM method introduces an additional term into the quasiparticle-phonon 
matrix elements, not taken into account in Ref. |p. In the present article we use 
the EOM form of the coupling part of the hamiltonian. In addition of being mi- 
croscopically more justified, this form of the coupling hamiltonian yields results 
closer to the experimental data and thus improves the quantitative predictibility 
of the MQPM. 

In the MQPM calculation we use for the protons the IpOf and 2sld0g oscillator 
major shells supplemented by the hii/2 intruder orbital from the next oscillator 
major shell. For the neutrons we use the 2sld0g7/2 and 2plf0h valence space. The 
proton single-particle energies 6^ of Eq. (1) are taken to correspond to the ones 
of Fig. 1, i.e. we take the same relative spacing of the key orbitals as used in 
the shell model calculation. The single-particle energies of Fig. 1 are extracted 
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Table 1 

Proton single-particle energies used in the MQPM calculations. The single-particle 
energies of the rest of the proton valence space are given in Fig. 1. 

s.p. orbital ^^Sj 137^3 i39La ^^^Pr ^^^Pm ^'^^Eu 



0f7/2 


-12.289 


-12.114 


-11.945 


-11.782 


-11.619 


-11.469 


0f5/2 


-8.808 


-8.694 


-8.583 


-8.476 


-8.367 


-8.269 


lP3/2 


-7.597 


-7.443 


-7.294 


-7.150 


-7.005 


-6.872 


lPl/2 


-6.212 


-6.083 


-5.959 


-5.838 


-5.717 


-5.618 


0g9/2 


-5.402 


-5.308 


-5.217 


-5.131 


-5.043 


-4.967 


0g7/2 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 



Table 2 

Neutron single-particle energies used in the MQPM calculations. 

s.p. orbital ^^^I ^^'^Cs ^^^La ^^^Pr ^^^pj^^ 145^^ 



1^5/2 
0g7/2 
2Si/2 
ld3/2 

Ohii/2 

lf7/2 
2P3/2 

Ohg/a 
2pi/2 

lf5/2 



-7.192 
-7.156 
-5.402 
-5.151 
-4.226 
0.000 
1.358 
1.724 
2.166 
2.715 



-7.251 
-7.283 
-5.425 
-5.203 
-4.341 
0.000 
1.499 
1.602 
2.371 
2.792 



-7.302 
-7.399 
-5.443 
-5.249 
-4.444 
0.000 
1.624 
1.488 
2.549 
2.853 



-7.345 
-7.504 
-5.456 
-5.289 
-4.537 
0.000 
1.735 
1.379 
2.704 
2.903 



-7.382 
-7.600 
-5.466 
-5.325 
-4.621 
0.000 
1.834 
1.275 
2.841 
2.945 



-7.413 
-7.688 
-5.473 
-5.356 
-4.696 
0.000 
1.923 
1.176 
2.963 
2.979 



from the experimental Sb spectrum |T6|, with exception of the 2si/2 single- 



particle energy which has not yet been measured. We have used the same value 
for the 2si/2 single-particle energy as in the work of Sagawa et al. The 
proton energies outside this set of states, as well as the neutron single-particle 
energies, we take from the Coulomb-corrected Woods-Saxon potential with the 
parametrization of |T^. These single-particle energies are displayed in Tables |l] 
and 0. 
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MeV 



3 - 



2 - 



1 - 



- 



1/2+ 2.990 
11/2 - 2.792 
3/2+ 2.708 



5/2+ 0.962 



7/2 - 



0.000 



Fig. 1. Adopted single-particle energies for the proton orbitals 2si/2) 1*^3/2) l'^5/25 0^7/2 
and /111/2 ill MQPM and the shell model calculation. 

In the next step a correlated ground state and the excited states of the even-even 
reference nucleus are constructed by use of the QRPA. In the QRPA the creation 
operator for an excited state (QRPA phonon) has the form 



(4) 



where the quasiparticle pair creation and annihilation operators are defined as 



A^{aa';JM) 



A(aa'; JM) = 



JM 



and aaa' = vT+Sa 



Here the greek indices u denote phonon spin J and parity tt. Furthermore, they 
contain an additional quantum number k enumerating the different QRPA roots 
for the same angular momentum and parity. Thus u = {JuijT^ui, k^j}- 

For each value of the angular momentum and parity the spectrum of the even- 
even nucleus is constructed by diagonalizing the QRPA matrix containing the 
usual submatrices A (quasiparticle-quasiparticle interaction) and B (induced by 
correlations of the ground state) |]13|. The basis states in our quasiparticle-phonon 
calculation are constructed from the previously determined BCS quasiparticles 
and the QRPA phonons. In the MQPM we make the following ansatz for the 
states in an even-odd nucleus 



; jm) 



n 



+ 



E 

nj'a 



nj 



(5) 
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where |— ) denotes the QRPA vacuum of the even-even reference nucleus. The 
overlap matrix elements between the quasiparticle-phonon states and the matrix 
elements of the quasiparticle-phonon hamiltonian in this basis have the following 
form 



1 



t 



-^^4'E{^ ^ t] ^ (6) 

X {hVta + En + hVta' + En' - 2Ea) X^n'X^nf^an an' ■, 

where Vta denote the QRPA-phonon energies, and 
A {ana'n']j) = JaJa' X! 



3n' 3a 



'yoL VOL' 33a \^ci' 

• • J I an' an r, ^ an an' 

3n 3 Jc' I "2 



3' 



(^an ^an' ' (7) 



Here = X^„, - The same definition holds for Y. 

The interaction matrix elements between the one-quasiparticle and quasiparticle- 
phonon states have the following form: 

Qaflnj . iia\'\-) = Hpp{aa'nn'Ja) (uaUa'X^^, - VaVa'Y"^,) all, 

3n' a<a' 

-li^Yl Hhh{aa'nn'Ja) {vaVa'X^^, - UaUa'Y^^,) a J 

3n' a<a' 

+ li^ J2 Hpi,{aa'nn'Ja) (uaVa'X^^, + VaUa'Y^^,) (8) 

3n' a<a' 

-li^Y. Hhpiaa'nn'Ja) {vaUa'X^^, + UaVa'Y^^,) 

3n' a<a' 



where 



Hpp{nn'aa'J) = 2vnUn'G{nn' aa' J) , (9) 

H-hhijin aa J) = 2unVn'G{nn' aa J) , (10) 

Hph{nn'aaJ) = 2vnVn' F{nn'aa' J) + 2u„M„/F(n'naa' J) (-!)■'" , (11) 

HY,p{nn'aa'J) = 2unUn'F{nn'aa' J) + 2t;„t;„/F(n'naa'J)(-l)^"+^"'+'^. (12) 

In the previous version of our model |^ the second term in the quasiparticle- 
phonon matrix elements of the hamiltonian in Eq. (6) was missing. This additional 
term stems from the use of the equations-of-motion (EOM) method |T5| when 
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deriving the eigenvalue equation (13) below. It has an important effect on the 
location of the three-quasiparticle-type states relative to the one-quasiparticle- 
type ones, and it is essential for yielding theoretical results in agreement with 
data. 

The overlap between the one-quasiparticle and the quasiparticle-phonon states 
is always zero. However, the overlap between two quasiparticle-phonon states 
can be non-zero and the quasiparticle-phonon states form a non-orthogonal over- 
complete basis set. The ansatz (5) leads to a generalized hermitian (or real and 
symmetric) eigenvalue problem which has the form |^ 

X! ^ij^f^ = -^n X! SijCf'\ (13) 
i 3 

where Hij = {i\H\j) and Sij = is the overlap matrix element between 

two basis states (one-quasiparticle or quasiparticle-phonon states). To solve this 
rather involved eigenvalue problem we adopt the method where we first solve the 
eigenvalue equation for the overlap matrix S: 



3 

The eigenvectors can be written in the basis as 



5: S.,«f = n,«f \ (14) 



They have the property of being mutually orthogonal, have a norm equal to 
unity and form a complete set after removing states having eigenvalue = 
(this removes the overcompleteness of the set {|«)}). 

Using the new orthogonal complete set of states (15) we can transform (13) to 
an ordinary real and symmetric eigenvalue problem of the form 



j:(i\H\j)9r = ^ngr\ m 

3 

where 

m~3) = ^^T.^f* {k\H\l)n? . (17) 

The coefficients of the eigenstates are calculated from the g coefficients in the 
following way: 

C: = Y.rik"'9t^uf\ (18) 

k 

In practice one omits states having eigenvalue less than some set upper limit 
e. 
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In the present work we have constructed the odd-A isotones by adding a proton 
particle to the adjacent even = 82 isotone. The even = 82 isotones thus 
provide our quasiparticles and QRPA phonons to be used in the ansatz of Eq. (5). 
As a matter of fact, these phonons have aheady been constructed in our previous 
study of the A^ = 82 isotones [Q. In the diagonalization of the MQPM matrix 
of Eq. (13) we have used the 4-6 lowest QRPA phonons of multipolarity 2+, 3", 
4"*", 5~, 6^ and 7~. This number of selected multipolarities is large enough to 
stabilize the spectrum of the odd isotones (see Ref. and is referring to the 
determination of the most relevant part of the three-quasiparticle hilbert space 
mentioned in the introduction. 

In this context it is important to point out that no 0"*" phonons are used to 
produce the MQPM results discussed in section ||, although in the first part 
of our work P the excited 0"*" states were discussed extensively. There are two 
reasons why we omit the O"*" phonons in the present MQPM calculation, namely, 
for the first, in the isotones ^^^Xe and ^^^Ba the calculated energies of the first 
excited O"*" state, 02", is seen to be more than 1 MeV too low in comparison 
with the corresponding experimental energy. In fact, in the QRPA the 0^ energy 
is lower than the 2f energy for these isotones. Thus the wave function of the 
0^ might not be adequately described and this particular phonon would produce 
two "spurious" low-energy excitations with spins 7/2"^ and 5/2"*" coming from the 
0^ ® 7/2"*" and 0^ ® 5/2+ quasiparticle-phonon couplings. Due to this unrealistic 
feature the 0"*" phonons are discarded in these isotones. Secondly, in the other 
isotopes under discussion the calculated energy of the state is well above the 
energy of the 2f state (confirmed by the experimental data in the cases of ^^°Ce, 
^^^Nd and ^"^^Sm) and calculations show that in this case the inclusion of the O"*" 
phonons has only a negligible effect on the results presented in section ^. 



3 Effective interaction and the shell model 



Here we will briefly sketch the theory of the effective interaction and the shell 
model (SM). 

Our scheme to derive an effective interaction can be divided into three steps. 
First one needs a free NN interaction V which is suitable for nuclear physics at 
low and intermediate energies. At present, the most viable approach seems to 
be the meson-exchange picture. Among the meson-exchange models, one of the 
more successful ones is the one-boson-exchange model of the Bonn group [[19 |. 



As the next step a reaction matrix G is introduced. In this way we overcome 
the problem that the NN potential has a strongly repulsive core which makes it 
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unsuitable for perturbative approaches. In this work we have calculated the G- 
matrix using the so-called double-partitioning scheme, see Ref. [10| for a recent 
review. 



The last step is to define a two-body interaction in terms of the G-matrix. We 
include all diagrams to third order in perturbation theory and sum the so-called 
folded diagrams to infinite order, an approach known as the folded diagram 
method, see e.g., Refs. PJT(][] for further details. 



Our basic approach in solving the many-body eigenvalue problem is the Lanczos 
algorithm. This is an iterative method which was first applied to nuclear physics 
problems by Whitehead |^ . The eigenstates are expanded in an m-scheme slater 
determinant basis, which implies that the dimension of the problem grows rapidly 
with increasing number of valence particles, see Table ^. The advantage of this 
representation is however the very efficient implementation of the computer code. 
This is a shell model calculation where no truncations of configurations are made. 



For more details of the shell model algorithm, see Ref. pT |. 
Table 3 

Number of basis states for the shell model calculation of the N = 82 isotones, with the 
1^5/2, 0(77/2, l'^3/2) 2si/2 and 0/iii/2 single particle orbitals. 



System 


Dimension 


System 


Dimension 


System 


Dimension 


i34Te 


36 


i39La 


108 297 


i^^Sm 


6 210 638 


135j 


245 


i40ce 


323 682 




9 397 335 


i36xe 


1 504 


141pj. 


828 422 




12 655 280 


137Cs 


7 451 


142Nd 


1 853 256 




15 064 787 


i38Ba 


31 124 


143pj^ 


3 609 550 


148Dy 


16 010 204 



The model space for the shell model calculation and the effective interaction is 
defined by the = 4 oscillator sheU (1^5/2, 0(77/2, 1(^3/2, 2si/2). In addition we 
have included the intruder 0/iii/2 orbital from the = 5 oscillator shell. Our 
model space consists of the proton orbitals outside the ^^^Sn core, ranging from 
the closed Z = 50, = 82 core to the Z = N = 82 core. There are no neutron 
degrees of freedom involved in this model. These degrees of freedom are accounted 
for by the various terms of the perturbation expansion used to derive the effective 
interaction. The adopted single-particle spectrum is as displayed in Fig. |l|. 

Before we present our comparison between the shell model results and those 
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obtained with the MQPM, we would hke to draw the attention to differences 
between the two methods. From the discussion in section ^ we recall that the 
MQPM method employs exactly the same G-matrix interaction as the one used 
to derive the shell model effective interaction. The single-particle energies for pro- 
tons in the orbitals lc?5/20(77/2l(i3/22si/20/iii/2 defining the shell model space are 
the same as those used in the MQPM. In addition, no phenomenological adjust- 
ments are made of the G-matrix in the MQPM approach. However, the reader 
should note that the single-particle basis for protons is larger for the MQPM, al- 
lowing thereby for proton core excitations across the Z = 50 shell gap. Secondly, 
also neutrons are active, yielding neutron core excitations across the = 82 shell 
gap. This might become important for the description of some low-energy collec- 
tive excitations of the even N = 82 isotones. In the shell model approach these 
degrees of freedom are supposed to be accounted for by terms included in the 
perturbative expansion of the effective interaction. Substantial differences in the 
two approaches may therefore reveal that such low-energy collective excitations 
are not accounted for in the shell model approach. 



4 Results and discussion 

In Tables we present our results for the calculated energy spectra of the odd 
= 82 isotones of interest in this article. In these tables we also present the avail- 
able experimental data in the two first columns of the tables. It is to be noted 
that the spin assignments on the first column are experimental ones and the spins 
in parentheses are only tentative. The MQPM and the SM results are always pre- 
sented in the last two columns of the tables. The three middle colums are reserved 
for presenting the results concerning the extreme quasiparticle-phonon picture. In 
this picture the energy of the multiplet emerging from angular-momentum cou- 
pling of the n:th phonon of spin J and parity vr, J^, with the single-quasiparticle 
state j'" is obtained by simply summing the phonon and the quasiparticle ener- 
gies. The column "proposed config." lists the phonon-quasiparticle statesQ which 
are among the leading ones in the wave function of the MQPM. 

In the first subsection below we discuss the differences between the extreme quasi- 
particle picture and the MPQM. Thereafter, the MQPM results are compared to 
the corresponding results obtained from an unrestricted shell model calculation 
with the orbitals ld5/20gY/2ld3/22si/20hii/2 defining the shell model space, having 
in mind the discussion in section |[ 

^ In the extreme quasiparticle-phonon picture we assume that only the proton- 
quasiparticle states are active in the low-energy spectrum of the odd isotones. 
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4-1 Microscopic quasiparticle-phonon model 



In the extreme quasiparticle-phonon picture the two energies, {E{phen) and 
£^(calc) in Tables ^[|^, are obtained in the following way. For the phenomeno- 
logical energy, E{phen), we have taken the phonon energy from experiments (i.e. 
from the measured spectrum of the even-even N = 82 isotone with one less pro- 
ton) and the quasiparticle energy was deduced from the proton single-particle 
energies of Fig. |I] using the BCS expression 



where is the single-particle energy of the active proton orbital and Ap is the 
chemical potential for the protons. The chemical potential we have taken simply 
to correspond to the energy of the last occupied proton orbital in the even- 
even = 82 reference nuclei. One exception is the case of ^^^Pr where in the 
corresponding reference nucleus ^^°Ce the proton gj/2 orbital is completely filled 
and the chemical potential thus lies somewhere between the g7/2 and 6.^/2 orbitals. 
In this case one can deduce the value of the chemical potential by using the above 
formula to reproduce the experimental energy difference between the 7/2+ and 
5/2"'" (single-quasiparticle) states in ^"^^Pr. The quantity A in Eq. (p!9|) is the 
pairing gap and assumes the phenomenological value 



The calculated extreme quasiparticle-phonon energy, _E'(calc), was obtained by 
summing the QRPA-calculated energy (see Ref. |0) of the J'" phonons of the 
even reference nuclei with the single-quasiparticle energy coming from the BCS 
calculation. 

In the following we shall comment on the MQPM results shown in Tables The 
emphasis lies on analyzing the lowest-lying quasiparticle-phonon multiplets of the 
odd isotones with respect to their span in energy, their centroids and the relative 
location of the various spins within the perturbed multiplet. The multiplet of the 
extreme quasiparticle-phonon picture looses its energy degeneracy through the 
action of the residual interactions and the resulting multiplet we call a perturbed 
one. It is of interest to see how the breaking of the degeneracy evolves from 
nucleus to nucleus and what its characteristic features are. 

We start first by noticing that the residual interactions yield relatively little 
perturbation on the energies of the single-quasiparticle-type states (denoted by 
the symbol (sqp) in the tables) in the MQPM calculation. The agreement in the 
single-quasiparticle energies between the MQPM and experiment is good due to 
the success of the BCS calculation in describing the lowest excitations of the odd 




(19) 



A = 12A~^/2MeV. 



(20) 
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Table 4 

Low-lying states for ^^^I. Experimental angular momentum values in parentheses are 
tentative. Energies are given in MeV. The symbol sqp in the leading configuration 
column denotes a single-quasiparticle state. For further information concerning the 
extreme qp-phonon picture, see the text. 







Extreme qp-phonon picture 






U 


J2j \ KiX.\J J 


proposed config. 


£'(phen) 


£^(calc) 






7/2+ 


0.000 


7/2+fsqD) 


0.00 


0.00 


0.000 


0.000 


(5/2+) 


0.604 


5/2+ (sqp) 


0.38 


0.70 


0.649 


0.508 


(5/2+) 


0.870 


2+ 7/2+ 


1.28 


1.51 


1.405 


0.888 






2+ 7/2+ 


1.28 


1.51 


1.498(9/2+) 


1.349(9/2 






2+ 7/2+ 


1.28 


1.51 


1.533(3/2+) 


1.053(3/2 


11/2+ 


1.134 


2+ 7/2+ 


1.28 


1.51 


1.543 


1.322 






2+ 7/2+ 


1.28 


1.51 


1.560(7/2+) 


1.619(7/2 


15/2+ 


1.422 


4+ 7/2+ 


1.57 


1.85 


1.886 


1.683 


17/2+ 


1.994 


6^ 7/2+ 


1.69 


1.94 


2.033 


1.966 



isotones. This, in turn, means that the phenomenological mean field, represented 
by the proton single-particle energies of Fig. is consistent with the residual 
pairing interactions used in the present calculation. 

As a general feature of the MQPM results one can say that the MQPM describes 
better the heavier iV = 82 isotones (both even |]^ and odd) since the heavier 
isotones are not in the immediate vicinity of the Z = 50 closed core and thus the 
quasiparticle description is expected to be better justified. 

For ^^^I the wrong theoretical 2+ energy in ^^"^Te (see Ref. ^j) leads to too high an 
energy for the multiplet 2+®7/2+ (see Table ^) and thus to too high an energy for 
the associated MQPM states (i.e. to too high a centroid of the MQPM multiplet). 
The same is valid for the 4+ energy and the associated 4+ (8> 7/2+ multiplet. For 
the other even reference isotones the phonon energies of the 2+,3j~,4+ and 6+ 
states are in reasonable agreement with experiment and in the case of ^'^'^Ba, 
^^°Ce and ^'^^Nd the agreement is rather good 0]. 

Let us now discuss the lowest observed multiplets of the odd isotones, namely 
2t ® 7/2+ for A = 139 and 2^ ® 5/2+ for A = 141, 143, as well as 3r ® 5/2+ for 
^^^Eu. Since experimental data on the 2j'"(8>7/2+ multiplet is missing or incomplete 
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Table 5 

Low-lying states of ^^"^Cs. Legend as in Table ^ 







Extreme qp-phonon picture 








J"" 


£;(exp) 


proposed config. 


i?(phen) 


S(calc) 


£;(MQPM) 




7/2+ 


0.000 


7/2+ (sqp) 


0.00 


0.00 


0.000 


0.000 


5/2+ 


0.456 


5/2+ (sqp) 


0.38 


0.45 


0.404 


0.208 






2+ ® 7/2+ 


1.31 


1.50 


1.315(5/2+) 


0.994(5/2+) 






2+ 7/2+ 


1.31 


1.50 


1.471(9/2+) 


1.463(9/2+) 






2+ 7/2+ 


1.31 


1.50 


1.512(3/2+) 


1.233(3/2+) 






2+ 7/2+ 


1.31 


1.50 


1.529(7/2+) 


1.519(7/2+) 






2+ 7/2+ 


1.31 


1.50 


1.530(11/2+) 


1.422(11/2+) 


1/2+ 


1.490 


2+0 5/2+, 4+ ® 7/2+ 


1.69 


1.95-1.96 


1.923 


1.361 


9/2- 


1.868 


2+0 11/2-, 3r®7/2+ 


3.26-3.27 


3.10-3.62 


3.126 


2.061 


1/2+ 


2.150 


1/2+ (sqp) 


2.13 


2.23 


1.986 


2.093 



for ^'^^I and ^^^Cs, only ^^^La is left for comparison. From Table |^ one observes 
that the experimental and the MQPM centroids correspond to each other rather 
nicely and the width of both spectra is the same. In addition, the 3/2+, 5/2+, 7/2+ 
and 11/2+ members of the multiplet are reproduced by the MQPM rather well 
but the MQPM clearly fails for the 9/2+ state. In this case the phenomenological 
quasiparticle-phonon energy ii^(phen) = 1.44MeV is more or less the centroid of 
both the experimental and the MQPM multiplet. The corresponding calculated 
energy ii^(calc) = 1.50MeV is slightly higher. 

The above described features are also seen in the 2i ® 7/2+ multiplet in ^"^^Pr, 
where this multiplet is not any more the lowest one but above the 2^ ® 5/2+ 
multiplet (see Table |^. A clear agreement of the MQPM with the data, especially 
in the case of the 3/2+ member of the multiplet, is evident. 

The same type of analysis can be performed for the 2^ 5/2+ multiplet which 
is the lowest one in ^'^^Pr and ^^^Pm (see Tables ^ and || ). In Table |^ it is 
seen that for ^^^Pr the centroid of the experimental multiplet is well reproduced 
by the MQPM as is also the case for ^^^Pm. A closer look at the spectrum of 
Mipj, j-gyggi^ig the 3/2+ and 9/2+ states are too far away from the 7/2+ state 
and the states are more homogeneously distributed in the theoretical spectrum. 
The undisturbed quasiparticle-phonon energies, i?(phen) and -E'(calc), lie near 
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Table 6 

Low-lying states of ^^^La. Legend as in Table ^. The underlined spin assignments are 
considered to be favoured by the SM and the MQPM calculations. 









Extreme qp-phonon picture 










£'(exp) 


proposed config. 




S(calc) 


£;(MQPM) 




7/2+ 


0.000 


7/2+ (sqp) 


0.00 


0.00 


0.000 


0.000 


5/2+ 


0.166 


5/2+ (sqp) 


0.38 


0.19 


0.154 


0.061 


1/2+ 


1.209 


2+0 5/2+ (4+0 7/2+) 


1.82 (1.90) 


1.69 (2.08) 


1.633 


1.331 


9/2+ 


1.219 




5 7/2+ 


1.44 


1.50 


1.451 


1.432 


(5/2)+ 


1.257 


2^g 


5 7/2+ 


1.44 


1.50 


1.256 


0.983 


(9/2+) 


1.381 


2|g 


5 5/2+ 


1.82 


1.69 


1.719 


1.491 


(11/2)- 


1.420 


11/2 


-(sqp) 


1.95 


1.71 


1.660 


1.770 


5/2+, 7/2+ 


1.421 


2+g 


5 5/2+ 


1.82 


1.69 


1.644 


1.467 


(9/2+) 


1.476 


4+g 


5 7/2+ 


1.90 


2.08 


2.068 


1.818 


11/2+ 


1.534 


2+g 


5 7/2+ 


1.44 


1.50 


1.530 


1.389 


7/2+ 


1.538 


2+g 


5 7/2+ 


1.44 


1.50 


1.462 


1.530 


3/2+, 5/2+ 


1.558 


2^g 


5 7/2+ 


1.44 


1.50 


1.495 


1.037 


(15/2)- 


3.375 


2^ ® 


11/2- 


3.39 


3.21 


3.242 


3.193 


17/2- 


3.927 


4j^ ® 


11/2- 


3.85 


3.79 


3.311 


3.798 



the top of the experimental and theoretical multiplet, respectively. This situation 
is reversed in the case of ^^^Pm, namely now the experimental multiplet is more 
homogeneously distributed due to the fact that the MQPM fails in predicting 
the location of the 9/2+ state whereas for the other members of the multiplet 
the MQPM energies roughly correspond to the experimental ones even though 
the 5/2+ and 7/2+ states are inverted in the theory. Once again the undisturbed 
quasiparticle-phonon energies are near the top of the perturbed multiplet. 

Finally, in Table ^one can observe the lowest multiplet, 3j~ 05/2+, of ^^^Eu. One 
can see that the centroid of the calculated multiplet is more or less correct but that 
the multiplet is far too compressed when compared with the experimental span 
of the multiplet (in experiment the 9/2" state comes very much down in energy). 
It seems that in the MQPM the biggest qualitative and quantitative problems 
appear for the 3/2" and 9/2" members of the multiplet. It has to be noted that 
in Table ^ the labelling of the higher excited states by the extreme quasiparticle- 
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Table 7 

Low- lying states of ^^^Pr. Legend as in Table ^. 



Extreme qp-phonon picture 





£^(exp) 


proposed config. 


7—1/1 \ 

i^(pnenj 


£/(caic) 


7 — 1 / TV T ^ \ T TV !V\ 

i?(MQPM) 


-E(SM) 


5/2+ 


0.000 


5/2+ (sqp) 


0.00 


0.00 


0.000 


0.000 


7/2+ 


0.145 


7/2+ (sqp) 


0.145 1) 


0.072 


0.104 


0.008 


11/2- 


1.118 


11/2- (sqp) 


1.31 


1.31 


1.301 


1.497 


3/2+ 


1.127 


2+ 


5/2+ 


1.60 


1.48 


1.314 


0.936 


(5/2)+ 


1.293 


2+ 


5/2+ 


1.60 


1.48 


1.371 


0.993 


1/2+ 


1.299 


1/2+ (sqp) 


1.49 


1.42 


1.396 


1.148 


3/2+ 


1.436 


3/2+ (sqp) 


1.23 


1.30 


1.342 


1.140 


(7/2)+ 


1.452 


2^ 


5/2+ 


1.60 


1.48 


1.440 


1.452 


9/2+ 


1.457 




5/2+ 


1.60 


1.48 


1.522 


1.413 


11/2+ 


1.494 




7/2+ 


1.74 


1.55 


1.622 


1.418 


9/2+ 


1.521 




7/2+ 


1.74 


1.55 


1.549 


1.413 


(3/2)+ 


1.608 




7/2+ 


1.74 


1.55 


1.626 


1.372 


1/2+ 


1.658 




5/2+ 


1.60 


1.48 


1.602 


1.545 


13/2+ 


1.768 




5/2+ 


2.08 


2.14 


2.206 


1.939 


15/2+ 


1.797 




5/2+ 


2.11 


2.21 


2.278 


1.978 


13/2+ 


1.986 


6^ 


5/2+ 


2.11 


2.21 


2.206 


2.004 


17/2+ 


2.070 


6| 


5/2+ 


2.11 


2.21 


2.280 


2.064 



1) 



The chemical potential of Eq. (19) was fixed by the experimental quasiparticle energy. 

phonon wave function is to be taken with a grain of salt since many of the MQPM 
states at these energies can be constructed by superposition of a large number of 
possible quasiparticle-phonon components lying rougly at the same energy region. 
Thus the overlap of the proposed quasiparticle-phonon structures of Table || with 
the MQPM wavefunctions is not necessarily very big. The assignments may be 
considered more as a means of keeping track of the experimental states with the 
same spin. 

As it was stated earlier, the spreading of the multiplet of the extreme quasipar- 
ticle-phonon picture can be ascribed to the influence of the residual interactions, 
coming mainly from the i^ai part of the residual hamiltonian. This spreading in 
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Table 8 

Low-lying states of ^^'^Pm. Legend as in Table ^ The underlined spin assignments are 
considered to be favoured by the SM and the MQPM calculations. 



Extreme qp-phonon picture 





^(exp) 


proposed config. 


-E'(phen) 


£;(calc) 


£;(MQPM) 


E{SM) 


5/2+ 


0.000 


5/2+ 


(sqp) 


0.00 


0.00 


0.000 


0.000 


7/2+ 


0.272 


7/2+ 


(sqp) 


0.39 


0.32 


0.342 


0.003 


11/2- 


0.960 


11/2" 


(sqp) 


1.08 


1.01 


1.009 


1.203 


3/2+ 


1.060 


3/2+ (sqp) 


1.01 


1.06 


1.088 


0.728 


1/2+ 


1.173 


1/2+ (sqp) 


1.26 


1.14 


1.166 


0.731 


(3/2,5/2) 


1.287 


2^ 


5/2+ 


1.58 


1.46 


1.176 


1.391 


3/2+ 


1.403 


2+ 


7/2+ 


1.96 


1.79 


1.854 




9/2+ 


1.456 


2+ ® 


5/2+ 


1.58 


1.46 


1.531 


1.512 


3/2+, 5/2+ 


1.515 


2+ 


5/2+ 


1.58 


1.46 


1.524 


1.325 


(5/2+) 


1.566 


2+ (E) 


5/2+ 


1.58 


1.46 


1.521(7/2+) 


1.430(7/2+) 


(9/2)+ 


1.566 


2+ 


7/2+ 


1.96 


1.79 


1.747 




3/2+, 5/2+ 


1.614 


2^ (E) 


7/2+ 


1.96 


1.79 


1.680 


1.403 


11/2+ 


1.664 


2t 


7/2+ 


1.96 


1.79 


1.854 




1/2+ 


1.753 


2+ 


5/2+ 


1.58 


1.46 


1.537 


1.632 


15/2+ 


1.898 




7/2+ 


2.49 


2.46 


2.528 




17/2+ 


2.288 


6t ^ 


5/2+ 


2.21 


2.36 


2.432 





energy ranges roughly from 100 keV to 400 keV in the discussed odd isotones and 
the calculated extreme quasiparticle-phonon energy {E{calc) in the tables) lies 
rather close to the top of the perturbed multiplet. This means that the residual 
interactions tend to redistribute the multiplet energies towards lower energies 
(the centroid of the perturbed multiplet is always substantially lower than the 
extreme quasiparticle-phonon energy). 

It is to be noted also that in the MQPM calculation the energies of the 3/2+ 
and 9/2+ members of the 2f ^ 7/2~^ multiplet are the ones least affected by the 
residual interactions in the case of ^^^I, ^'^^Cs, ^^^La and ^^^Pr. For the 2f (8)5/2+ 
multiplet all the spin members of the multiplet are clearly affected by the residual 
interaction, some of them considerably. The most strongly affected member of the 
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Table 9 

Low-lying states of ^^^Eu. Legend as in Table |. In this case only the MQPM results 
are available since the SM calculations are already beyond any reasonable effort. 



£;(exp) 



Extreme qp-phonon picture 
proposed configuration(s) £'(phen) £'(calc) ii^(MQPM) 



5/2+ 
7/2+ 
11/2- 
1/2+ 
3/2+ 
9/2- 
1/2+ 
7/2~ 

3/2", 5/2" 
3/2- 
11/2- 



7/2- 
3/2+ 
3/2- 
5/2" 
11/2" 
9/2- 



1/2+, 3/2+ 
(3/2,5/2)+ 

3/2+ 

5/2+ 
9/2-, 11/2- 



0.000 
0.329 
0.716 
0.809 
1.042 
1.368 
1.460 
1.500 

1.567 
1.600 
1.602 



1.745 
1.758 
1.762 
1.766 
1.792 
1.827 



3/2+,5/2+ 1.845 



1.881 
1.915 
2.049 
2.114 
2.117 



5/2+ (sqp) 
7/2+ (sqp) 
11/2- (sqp) 
1/2+ (sqp) 
3/2+ (sqp) 
3f 5/2+ 
2+ 5/2+ 
3f 5/2+ 
3f 5/2+ 
3^ 5/2+ 
5/2+ 
5/2+ 



2l 



5/2^ 
2+ 5/2^ 
2+ 5/2^ 
3r®7/2+, 2t^ 
5/2^ 
3r 7/2^ 



11/2" 



2+ 
^1 



37 7/2+ 
3f 7/2+, 2+ 
3f 7/2+, 2+ 
2+0 7/2+, 3^ 
2+0 7/2+, 3^ 
2+ 7/2+, 3^ 



2^ 5/2^ 



11/2" 
11/2" 
11/2" 
11/2" 
11/2" 

11/2" 



2+0 7/2+, 3r 

2+ 7/2+ 
2+ 7/2^^ 

3f 7/2+, 2+ 11/2- 



0.00 

0.39 

1.08 

1.26 

1.01 

1.81 

1.66 

1.81 

1.81 

1.81 

1.81 

1.81 

1.66 

1.66 

1.66 
2.20,2.57 

1.66 

2.20 

2.20 
2.20,2.57 
2.20,2.57 
2.05,2.89 
2.05,2.89 
2.05,2.89 

2.42 
2.05,2.89 

2.05 

2.05 
2.20,2.57 



0.00 

0.54 

0.60 

0.76 

0.71 

1.48 

1.48 

1.48 

1.48 

1.48 

1.48 

1.48 

1.48 

1.48 

1.48 
2.02,2.08 

1.48 

2.02 

2.02 
2.02,2.08 
2.02,2.08 
2.02,2.07 
2.02,2.07 
2.02,2.07 

2.58 
2.02,2.07 

2.02 

2.02 
2.02,2.08 



0.000 
0.551 
0.618 

0.788 

0.737 

1.542 

1.547 

1.541 
1.532(1/2") 
1.542(5/2-) 

1.535 

1.553 
1.543(9/2+) 
1.571(7/2+) 
1.587(5/2+) 

2.080 

1.258 

2.070 

2.065 

2.084 

2.084 
1.928(7/2+) 
1.952(5/2+) 
1.977(9/2+) 
2.090(1/2+) 
2.161(5/2+) 

2.295 

2.252 
2.095(9/2-) 



2f (g> 7/2+ multiplet is the 5/2"*" state whereas 3/2"*" state is most affected in the 
2|®5/2+ multiplet. 



4-2 Comparison with shell model results 



In the MQPM scheme one starts with the BCS occupation probabilities stemming 
from the even isotones, given in Ref. [^]. As discussed in section 0, the residual 
many-body hamiltonian can be expressed in terms of these occupation probabil- 
ities. In general, for the single-particle states we obtained a good agreement in 
Ref. between the shell model occupation probabilities and those of a BCS cal- 
culation for the even isotones. Since the states of the MQPM, discussed in Tables 
^"P, are combinations of single-particle states and phonons from the correspond- 
ing even isotones, it is then of interest to see whether one can retrace eventual 
discrepancies between the shell model approach and the MQPM to the fact that 
the SM employs a smaller set of single-particle orbitals in the diagonalization 
than the MQPM. It is therefore important to see how a renormalized effective 
interaction obtained by perturbative many-body methods is able to account for 
degrees of freedom not accounted for by the SM model space. It is, however, im- 
portant to notice that although the MQPM employs a larger set of single-particle 
orbitals, only a limited set of states are obtained from the diagonalizations. In 
the SM all states are, in principle, taken into account in the diagonalization. 



The shell model results are presented in Tables In, general, there is a fairly 
good agreement between the shell model results and experiment, with deviations 
of the order of 100—300 keV. However, as the number of valence particles increases 
from 3 in ^^^I to 13 in ^^^Pm, the description of the spacing between the lowest- 
lying 5/2+ and 7/2+ states gets worse. 

In order to understand these differences we have performed additional shell 
model calculations for ^^^Pr with just the G- matrix as effective interaction, in 
order to see whether different approaches to the effective interaction within the 
2sld0gY/20hii/2 model space yield significant discrepancies. These results are dis- 



played in Table [1^, under the column labelled E{G). The shell model results 
with the effective interaction to third-order in G from Table 0, together with the 
corresponding MQPM results, are included for comparison. These are labelled by 
^(SM) and E(MQPM), respectively. 



One can note from Table |10| that when going from the third-order effective in- 
teraction over to the bare G-matrix interaction in the shell model, the spacing 
between the lowest-lying 5/2+ and 7/2+ states increases to 0.034 MeV, although 
it is still far from the experimental value. The spectrum of the other states is 
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more compressed than for the £'(SM) results. The proton occupation numbers 
do also change, but all states of ^^^Pr are still strongly dominated by admixtures 
from the 5/2"*" and 7/2"*" single-particle states. 

Table 10 

Lowest-lying states of given multipolarity in i^^Pr using various approximations to the 
effective interaction. 





^(exp) 


S(BCS-l) 


^(BCS-2) 


S(MQPM) 


E{G) 


S(SM) 


5/2+ 


0.000 


0.17 


0.00 


0.00 


0.00 


0.000 


7/2+ 


0.145 


0.00 


0.07 


0.104 


0.034 


0.008 


11/2" 


1.118 


1.42 


1.31 


1.301 


0.889 


1.497 


3/2+ 


1.127 


1.44 


1.30 


1.314 


0.915 


0.936 


1/2+ 


1.299 


1.58 


1.42 


1.396 


1.292 


1.148 



The question then arises whether these differences can be traced back to the 
use of a smaller model space in the shell model calculation. We have therefore 
performed a BCS calculation in the model-space defined by the 2sl(i0(77/20/;.ii/2 
single-particle orbitals, using the G-matrix as interaction and the single-particle 
energies of Fig. ^ These results are labelled i?(BCS-l) in Table |TU|. It has to 



be noted that in this BCS calculation the mean-field part is the same as in the 
E{G) calculation but that a BCS calculation includes only a very restricted set 
of states as compared to the full shell model diagonalization. The comparison 
reveals then how good or bad the BCS approximation is. As stated earlier, the 
BCS calculation forms the basis for the MQPM method. Moreover, the MQPM 
method used here employs a larger single-particle basis than that used in the 
shell model calculation or the BCS-1 calculation. It may therefore be of interest 
to see how the BCS calculation changes when we go to the model space employed 
in the MQPM calculation. The results of such a BCS calculation are denoted by 
E(BCS-2) in Table 0. 

One can see from Table that the results of the BCS-1 calculation clearly devi- 
ate from the BCS-2 results, indicating the importance of a larger single-particle 
basis. The relative positions of the 5/2+ and 7/2+ states are inverted and the 
other states are higher up in excitation energy in the BCS-1 calculation. Com- 
paring the BCS-1 results with the E{G) results, which are actually in reasonable 
agreement with experiment, clearly indicates that the simple BCS picture is far 
from sufficient in the restricted 2sl(i05f7/20/iii/2 valence space. The BCS-2 calcula- 
tion, which employs the larger set of single-particle states, yields results which are 
closer to experiment and close to the MQPM results. The mere difference between 
the BCS-2 and MQPM results is in the better reproduction of the 5/2+ — 7/2+ 
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spacing by the MQPM. 



In summary, the results in Table |10| seem to indicate that degrees of freedom 
not accounted for by the model space employed in the shell model calculation, 
are important in order to get a proper reproduction of the experimental spacing 
between the two lowest-lying 5/2+ and 7/2+ states. However, when comparing 
the BCS-1 results with those of a shell model calculation with the bare G- matrix, 
one sees that there are important differences. Sources of these discrepancies are 
the many-body configurations not accounted for by the BCS approach. How these 
differences will appear in a shell model calculation which would employ the same 
set of single-particle energies as the BCS-2 or MQPM approaches, is however not 
clear. 

In column three of Tables]^ - ^ we indicate the MQPM states that are supposed to 
be of single-quasiparticle nature. For comparison with the shell model we create 
single-quasiparticle states for the odd-nucleon system by coupling a (j'^) particle 
to the 0+ ground state of the neighbouring A — 1 even system, a]\SM{A - 1)). 
By calculating the squared overlap between the constructed single-quasiparticle 
state and the SM state, Kj'^; SM(A)|a]|0+; SM(A - we obtain a measure 

of the fraction of single-quasiparticle structure in our SM wave function. The 



results are tabulated in Table |TT|. There is a nice correspondence between the 
Table 11 

Squared overlaps \{r;SM{A)\a]\0-^]SM{A - for the SM and one-quasiparticle 
probabilities for the MQPM. 

jn 1351 137cs 139La ^^^Pr ^^^p^^ 

SM MQPM SM MQPM SM MQPM SM MQPM SM MQPM 



7/2+ 


0.99 


0.99 


0.99 


0.99 


0.98 


0.99 


0.97 


0.99 


0.95 


0.99 


5/2+ 


0.94 


0.98 


0.95 


0.98 


0.95 


0.99 


0.94 


0.98 


0.98 


0.98 


3/2+ 


0.00 


0.01 


0.01 


0.03 


0.07 


0.19 


0.12 


0.45 


0.89 


0.97 


3/2+ 


0.00 


0.00 


0.09 


0.00 


0.02 


0.00 


0.75 


0.52 


0.00 


0.01 


1/2+ 


0.11 


0.00 


0.11 


0.15 


0.29 


0.32 


0.89 


0.73 


0.90 


0.97 


ii/2r 


0.98 


0.94 


0.96 


0.95 


0.94 


0.95 


0.91 


0.96 


0.89 


0.96 



MQPM states proposed to be of single-quasiparticle nature and those SM states 
with predominantly single-quasiparticle structure. In the case of the 3/25*" and 
3/2^ states in ^^^Pr there is a strong mixing of the one- and three-quasiparticle 
components in the MQPM, stronger than in the SM. Squared overlaps of mag- 
nitudes 0.95 - 1.00 confirm that the 5/2+ and the 7/2i states are fairly pure 
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single-quasiparticle states. 



We then turn to a comparison of the predominantly three-quasiparticle states in 
the two models. For this we consider the multiplets 2^07/2"'' for A = 135, 137, 139 
and 2t 5/2+ for A = 141, 143. We start with the 2^ (g) 7/2+ multiplet. Looking 
at Tables one can see differences in the SM and the MQPM spectra of states 
(3/2+, 5/2+, 7/2+, 9/2+, 11/2+) belonging to this multiplet. Characteristic fea- 
tures are 1) the centroid of the SM multiplet is always lower than the centroid of 
the MQPM multiplet. 2) The high-spin members of the multiplet (7/2+, 9/2+, 
11/2+) correspond quite well to each other in the SM and the MQPM whereas 
the 5/2+ member of the multiplet is consistently lower in the SM spectra than 
in the MQPM ones (in both cases the 5/2+ state is usually the lowest one). 3) 
Differences show up for the lowest-spin member of the multiplet, i.e. for spin 
3/2+, which is always the second lowest in energy in the SM spectra but among 
the three highest levels in the MQPM spectra. 

The above listed properties of the SM and the MQPM spectra distinguish between 
the two calculations. For ^'^^I and ^^''Cs the data is missing or incomplete so that 
essentially only ^^^La is left for comparison. From Table |^ one observes that the 
experimental and the MQPM centroids correspond to each other rather nicely 
and the width of both spectra is the same. In addition, the 3/2+, 5/2+, 7/2+ and 
11/2+ members of the multiplet are reproduced by the MQPM rather well (both 
the SM and the MQPM fail for the 9/2+ state) whereas in the case of the SM 
large deviations are observed for the 5/2+ member and especially for the 3/2+ 
member of the multiplet. Thus the experiment favours the MQPM sequence of 
levels in the 2i ® 7/2+ multiplet (although also in the MQPM the sequence of 
levels is not completely correct). For the 2f ® 7/2+ multiplet in ^^^Pr (see Table 
1^) the differences between the two models and the experimental data are small. 

The same type of analysis can be performed for the 2]*" eg) 5/2+ multiplet which is 
the lowest one in ^'^^Pr and ^^^Pm (see Tables |^ and ^ . In this case the differences 
between the SM and the MQPM spectra are less than for the 2+ ® 7/2+ case. 
For Pr the centroid of the experimental multiplet is well reproduced by the 
MQPM, clearly better than by the SM, whereas for ^^^Pm both the SM and the 
MQPM have roughly the correct centroid. The largest difference between the 
SM and the MQPM spectra is found in the location of the 5/2+ member of the 
multiplet and the experimental data favours the MQPM for the 5/2+ energy. 
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5 Conclusions 



The present work discusses tlie tlieoretical interpretation of low-energy excita- 
tions of odd = 82 isotones between the mass numbers A = 135 and A = 143. 
The energy spectra of these isotones have been calculated by using the micro- 
scopic quasiparticle-phonon model (MQPM) and the results have been compared 
with the extreme quasiparticle-phonon picture and the results of a large-basis 
shell model calculation with 3 — 11 valence protons outside the doubly- magic 
^^^Sn core. This work is a direct continuation of our earlier work on even N = 82 
isotones and the same realistic, microscopic two-body G-matrix interaction 
has been used in the present MQPM calculation as was used in the QRPA calcula- 
tion of the even isotones (in fact, the QRPA calculation is a necessary prerequisite 
of the present MQPM calculation). Also in the shell model calculation we use the 
same effective two-body matrix elements, derived from the above-mentioned G- 
matrix elements through many-body perturbation techniques, which were used 
for the even isotones. 

Overall, the spectra of the odd isotones are described well by the MQPM consider- 
ing that no fitting of the interaction was done. This feature may be traced back to 
the capability of the BCS approach in describing excitations of one-quasiparticle 
type. From this one can conclude that the pairing-type of residual interactions 
used in the MQPM are consistent with the mean field extracted from the experi- 
mental single-particle energies in ^^^Sb. The low-lying experimental levels can be 
labeled by their assumed leading quasiparticle-phonon contributions and their en- 
ergies can be compared with the energies of the unperturbed quasiparticle-phonon 
multiplets and the perturbed ones emerging from the MQPM calculations. In the 
MQPM both the widths and the centroids of the perturbed multiplets, as well 
as the sequences of different spins within the multiplets, are described rather 
nicely. The shell model describes most individual states very well, in many cases 
better than the MQPM, but has difficulties in describing the centroids of the 
quasiparticle-phonon multiplets and energies of some members of these multi- 
plets, particularly the 3/2+ states in the 2f ® 7/2+ multiplets and the 5/2+ 
states in the 2+0 5/2+ multiplets. 

The aim of the present work was to see how well a truly microscopic model, based 
on quasiparticle-phonon coupling and realistic microscopic G-matrix interactions, 
can describe the level systematics of a set of heavy semi-magic nuclei. At the same 
time the results of these calculations can be compared with results coming from 
a large-scale shell model calculation with truly microscopic effective interaction 
based on the same G-matrix which is used in the quasiparticle-phonon calcula- 
tion. Considering that in both calculations only very few parameters enter the 
calculation, the success of both models is surprisingly good. 
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